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Abstract 

In this paper we study Peierls instabilities for a half-filled two-dimensional 
tight-binding model with nearest-neighbour hopping t and next nearest- 
neighbour hopping t' at zero and finite temperatures. Two dimerization pat- 
terns corresponding to the same phonon vector (ir, tt) are considered to be 
realizations of Peierls states. The effect of imperfect nesting introduced by t' 
on the Peierls instability, the properties of the dimerized ground state, as well 
as the competition between two dimerized states for each t! and temperature 
T, are investigated. It is found: (i). The Peierls instability will be frustrated 
by t' for each of the dimerized states. The Peierls transition itself, as well 
as its suppression by t! ', may be of second- or first-order, (ii). When the 
two dimerized states are considered jointly, one of them will dominate the 
other depending on parameters t' and T. Two successive Peierls transitions, 
that is, the system passing from the uniform state to one dimerized state and 
then to the other take place with decrease of temperature for some t' values. 

Implications of our results to real materials are discussed. 
PACS numbers: 71.45.Lr, 63.20.Kr 
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I. INTRODUCTION 



The Peierls instability, associated with bond-order or charge- density waves, is an impor- 
tant phenomenon in low dimensional materials driven by electron-phonon interaction [|l|. 
The occurance of such an instability is directly relevant to the nesting property of the Fermi 
surface (FS) in the normal high temperature state of materials. It has been observed in quasi- 
one dimensional (ID) metals such as polyacetylene (CH) Z , blue bronzes A0.3M0O3 (A=K, 
Rb, Tl) ||, as well as quasi-two dimensional (2D) materials such as purple bronzes AM06O17 
(A=Na, K, Tl) and a series of monophosphate tungsten bronzes (P0 2 )4(W0 3 ) 2m with 
4 < m < 14 Hf]. 

Theoretically the description for the Peierls instability was mainly developed for ID 
metals, and relatively few investigations are concerned with 2D systems [fT0| ■ Actually 
the 2D study of the Peierls instability is not straightforward because of the more complex 
FS structure as well as richer distortion patterns. For example, in Fig. |]we have shown two 
possible dimerization patterns for the lattice distortion and the corresponding bond hopping 
in two dimensions. Both of them correspond to phonons with wave vector (tt, it), and their 
difference is that for pattern (a) the dimerization is in both directions, while it is only in 
one direction for pattern (b). The simplest model to discuss the Peierls instability in two 
dimensions is a square lattice tight-binding model with nearest-neighbour (n.n.) hopping 
t, i.e., the 2D version of the well-known Su-Schrieffer-Heeger (SSH) model |i~6| . It was 
originally studied by several authors in connection to high-T c superconductors [jT0|JTTH . At 
half-filling the FS is perfectly nested with nesting vector Q = (it, 7r) so that arbitrarily small 
electron-lattice coupling [|l5j will induce a Peierls instability into any one of the dimerized 
states shown in Fig. p] (if quantum lattice fluctuations are ignored). 

The above situation, i.e., perfect nesting of FS is, however, special because it may be 
easily broken for example, by introduction of next nearest-neighbour (n.n.n.) hopping t' . 
Actually for real materials, e.g., quasi-2D high-T c superconductors n.n.n. hopping was found 
to be important. Also for those quasi-2D materials which show a Peierls instability, perfect 
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nesting of their Fermi surfaces is not satisfied. It is reasonable to adopt a tight-binding model 
with n.n. and n.n.n. hopping to simulate an essential property of them: the relevance of 
nesting. Then a question is naturally asked: how is the Peierls instability affected by £'? 

It is the purpose of this paper to make a systematical study of Peierls instabilities for 
a 2D t-t' model, where the nesting property of the FS can be modulated. In contrast, the 
corresponding problem in the ID case might be less attractive because the FS is always 
composed of two points and perfectly nested as long as t'/t < 0.5, and then at least in this 
region t' has no effect on Peierls instability |l7j. The 2D problem at zero temperature was 
addressed in an earlier work ||15|| . It was found that the Peierls instability will be suppressed 
with increasing t' for each of the two patterns (a) and (b), when they are studied individually. 
Moreover, when they are considered jointly, one of them was found to dominate the other 
or reverse when t' is less or greater than some critical value t' ab , i.e., there exists a crossover 
between the two patterns at t' = t' ab . Further one may conjecture that the same behavior 
should also occur at finite temperatures, and most probably the value t' ab will change with 
temperature. Then an intriguing phenomenon may arise: at some t', the two dimerized states 
will dominate successively with change of temperature so that double Peierls transitions (one 
is between uniform state and dimerized state, the other is between two dimerized states) 
are expected to take place as a function of temperature. Actually, two successive Peierls 
transitions were observed in TlMo 6 Oi 7 || and most of the series (P0 2 )4(W0 3 ) 2m So 
it is interesting to extend the previous work to finite temperatures to discover the above 
phenomenon. At the same time, new findings at ground state will be supplemented. 
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FIG. 1. The lattice distortion patterns (a) and (b). In the figure a thick solid line corresponds 
to a strong bond with hopping integral t(l + 5), a dashed line corresponds to a weak bond with 
hopping integral t(l — 5), and a thin solid line corresponds to a normal bond with hopping integral 
t. Both patterns correspond to phonons with wave vector (vr,7r). The dimerization is along two 
axes for pattern (a), while only along the x axis for pattern (b). 

Before entering the details we implement a remark on the lattice distortion. Recently, 
Ono and Hamano studied the 2D Peierls instability numerically in the absence of t' and found 
that a complex multi-mode distortion pattern has a lower ground state energy than any one 
of the (it, 7r) structure presented here As also described by the authors themselves, 

however, such a distortion pattern is not unique and may be infinitely degenerate. So 
it is difficult to extract an explicit pattern for real applications. Moreover, no further 
confirmation on these new patterns is present at the moment. Thus the question of the 
real lattice distortion for the 2D Peierls system is still open even in the case of t' — 0. Under 
these circumstances we prefer to adopt the two commonly assumed dimerization patterns 
(a) and (b) as candidates of the Peierls state. These two patterns comply with the nesting 
vector Q = (jr, it) and realize an unconditional Peierls instability at t' = 0. Similarly, 
one may also argue that for each t' and temperature there may exist a lattice distortion 
structure, though perhaps extremely complicated, with the lowest (free) energy. We have 
no intention to include all kinds of lattice distortions, which is actually impossible in an 
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analytical treatment on an infinite lattice. Instead, throughout the work we constrain the 
lattice distortions to the patterns (a) and (b) (although a discussion on a generalization to 
these two patterns is given in the last section), and based on this, we intend to investigate 
the presence of t' after fully understanding the t' — case. This is a natural topic by physical 
consideration as illustrated in the preceding paragraphs, which is helpful to understand the 
2D Peierls instability itself as well as relevant materials even if the assumed dimerized states 
(a) and (b) might be not the lowest energy Peierls states. 

II. FORMULATION 

Our Hamiltonian represents a half filled 2D t — t' model including lattice displacements. 
It reads @: 

H = ~t J^i 1 + a ( U tj ~ U i+l,j)]( C tj,aCi+l,j,a + h.C.) -tJ2i 1 + a ( U ij ~ U L+l)K C L> C ^+l>* + hX -) 

K 

E ( c h,*Ci+id+i,° + cl h vd+ij-i,o + h.c.) + — Y,'(<, - u i+i,j? + Kj - u h+if} ■ (!) 

All above notations are conventional: cJ Ji0 .(cjj jCT ) is the creation (annihilation) operator for 
an electron at site with spin a (i denotes x coordinate and j denotes y coordinate); 
Ui j V is the displacement component of site in x/y direction; t, t' are n.n and n.n.n. 

hopping parameters and a is the electron- lattice coupling constant. The last term above 
describes the lattice elastic potential energy with K the elastic constant. The phonons are 
treated in adiabatic approximation. 

The lattice distortion, as shown in Fig. [I], may be explicitly expressed as 

for pattern (a) and 

for pattern (b), where u is the amplitude of dimerization which needs to be determined. 
Experimentally these two patterns may be differentiated by high resolution x-ray or neu- 
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tron scattering. For convenience, two dimensionless parameters are defined as follows: the 
dimerization amplitude 5 = au and the electron-lattice coupling constant rj = a 2 t/K. 

The Hamiltonian (|1]) may be easily diagonalized to give the following electronic spectra 
in momentum space for patterns (a) and (b), respectively, 

£ t,a = ~ 4*' cos k x cos k y ± 2\J (cos k x + cos k y ) 2 + <5 2 (sin k x + sin k y ) 2 , 

e^ b = — 4t' cos k x cos ky ± 2\J (cos k x + cos A;,,) 2 + 5 2 sin 2 fc^ . (2) 

Before continuing we check the symmetries of the above spectra. It is easy to see that the 
spectra in both cases are invariant under the same operation: k x — > k x ±Tr and simultaneously 
k y —> Oky ± 7T. In addition, the spectra e^ a have reflection symmetry: k x — > —k x and 
simultaneously k y — > — as well as exchange symmetry: k x ^ k y ; while the have only 
the reflection symmetry: k x — > — fe^ or A;^ — > — k y , or both. 

The electronic grand canonical partition function S may be written as: 

lnH = 2^1n[l + e- /3(£ ^)] (3) 

with /3 = l/ksT (T: temperature) and ^ = ±: band index. The factor 2 in front of the 
summation includes spin degeneracy. The chemical potential fi is adjusted to yield the right 
filling, i.e., N e = ^JMnH. N e is the total electron number, equal to the total number of 
lattice sites iV at half-filling. 

The required optimal dimerization parameter 5* is determined by minimization of the 
total free energy which is given by 

F = -^InE + N ef i + E L . (4) 

Here El denotes the lattice elastic energy which is, in unit of t, N5 2 /r] for pattern (a) and 
NS 2 /(2r]) for pattern (b). Note that El for (a) is twice of that for (b) under the same 5. 

Throughout the paper we take the Bolzmann constant fee = 1 and the n.n. hopping 
integral t as the energy unit. The parameter 77 is fixed at a typical value 0.5. The variation 
of rj should not change the results qualitatively, as has been seen at ground state in the 
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earlier work [15] . Also notice that the Hamiltonian with —t' may be mapped onto that with 
t' at half-filling through the transformation: Cjj )(J — > (— l) l+J q j a , so we only consider t' > 
in the following. 



III. RESULTS AT T = 



The effects of t' on the Peierls instability were carefully studied at ground state in Ref. 



| 15|j , and the main results are summarized here. For each of the two patterns, it was found 
that t' tends to suppress the Peierls instability. But the details in both cases are different: 
for (a) the suppression is of first-order while for (b) it is of second-order. The critical values 
of t' for suppression in both cases are close but slightly different. More interestingly, a 
crossover between the two dimerized states themselves may take place with change of t' . 

After a brief presentation of the above results, our main attention in this section is paid 
to other useful properties in the ground state, which were not discovered previously. First 
let us have a look at the variation of the FS before and after dimerization. This helps 
thoroughly understand the impact of the FS structure on the Peierls instability. We take 
t' = 0.15 for example and plot the Fermi surfaces in Fig. |2| for both patterns (also cf. 
the plots of electronic spectra shown later). In this figure, the FS before dimerization is 
shown by the gray curve in the original Brillouin zone, while after dimerization the Fermi 
surfaces for the ' — ' and '+' bands are shown in the reduced Brillouin zone by the solid and 
dashed curves, respectively. Please note by the way the differences of the FS structures for 
both patterns which are controlled by the symmetries of the respective spectra addressed 
above. From Fig. |^ a common observation in both patterns is that after dimerization there 
are a few electrons resting in the '+' band, whose energy levels are lifted as seen from Eq. 
(0). That is, the energies for this small fraction of electrons are increased by dimerization, 
leading to a partial cancellation of the energy gain acquired by most of electrons dropping 
into the '— ' band. This situation is in contrast to that for t' = where the original FS is 
perfectly nested so that all electrons drop into the '— ' band with each energy level reduced 
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(at least not raised) by dimerization. Therefore, one may see that the imperfect nesting of 
the original FS induced by t' is unfavorable to gain electronic energy to the largest extend. 
And once such a gain is not sufficient to compensate the cost of the elastic energy, the Peierls 
instability will be prohibited. 




FIG. 2. The Fermi surfaces before dimerization (gray curve) and after dimerization (solid and 
dashed curves) for t' = 0.15 for both patterns. For pattern (a), the FS of the '+' band (dashed 
curve) is formed just by the two pockets around (n/2,— n/2) and (— vr/2, vr/2); while for pattern 
(b) it is formed by the four pockets. The ' — ' band for each case has the large FS (solid curve) with 
nearly the square shape, except for the cut corners. The square: jA^I + \k y \ = ir is the reduced 
Brillouin zone (also the FS in the case of t' = 8 = 0). 

Another observation we want to illustrate is that the dimerized state (a) or (b), due to the 
Peierls transition, is metallic or semi-metallic in two dimensions even for t' = rather than 
insulating as in one dimension. Actually there is no finite gap opened by the dimerization 
in any case. To more clearly see this point we have plotted the electronic spectra in the 
dimerized states for several t' values in Fig. as well as their corresponding density of 
states (DOS), i.e., 4 J2\t,v $( £ ~ £ V) i n Fig- The two bands are found to touch or overlap 
for each t' and for each pattern, as also seen from DOS where no gap appears in each case. 
We point out that the Peierls instability, when occuring in those real quasi-2D materials, is 
always associated with a metal-metal transition ||. 

At this stage, it is worthwhile to mention several points seen from Fig. i. The 
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two bands in the region (— tt, 0)-(0, it) are identical for pattern (a). (Thus they are always 
overlapped.) ii. The spectra along the path (— %, 0)-(0, 7r)-(7r, 0) are symmetric about point 
(0, 7i") for pattern (b), while it is not for pattern (a). This can be also seen from Fig. 0. hi. 
For pattern (b), the two bands display no overlap for t' — 0.1 (in contrast to the case for 
t' = 0.15). Then the ground state energy is simply a summation for the ' — ' band, which 
becomes t'- irrelevant. So the optimal parameter 5* should take the same value as that for 



t' = 0, as discussed in the earlier work [15 
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FIG. 3. The electronic spectra in the dimerized states for several t' values for both patterns. 
The solid and dotted curves in each panel represent energy bands '— ' and '+', respectively, and 
the dashed horizontal line shows the chemical potential. Note that the two bands in the region 
(— 7T, 0)-(0, 7r) are identical for pattern (a). 
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FIG. 4. The density of states (DOS) for each case corresponding to Fig. ||[ The dashed vertical 
line in each panel shows the chemical potential. 

At last, we point out that the DOS in each dimerized state for finite t' exhibits the 
feature of asymmetric double peaks. This could provide characteristic information on some 
measurable properties like optical absorption in real materials when available. 



IV. PEIERLS TRANSITIONS AT T > 

We discuss finite temperature transitions in this section. For clearness we first look at the 
results for the two cases of pattern (a) and (b) separately, and then consider the competition 
between them. The optimal dimerization parameters with change of temperature were 
calculated for various values of t' . The results for pattern (a) are shown in Fig. |5| and those 
for pattern (b) in Fig. [| A notable common property seen from both figures is that the 
Peierls transition temperature (denoted as and for pattern (a) and (b), respectively) 
decreases with increase of t' . This just means that the n.n.n. hopping t' is unfavorable to 
the Peierls instability to occur, which is expected and consistent with the results at ground 
state. It deserves to be pointed out that the lowering of the transition temperature T p with 
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increasing t' is consistent with the experimental findings in the bronzes (P02)4(WOs)2m, 
where T p was found to increase with m because of better nesting properties ||. Then we 
continue to look at the types of transitions in both cases. For pattern (b), the optimal 
dimerization parameter 5* goes to zero smoothly with increase of T for each t', i.e., the 
Peierls transition is always of second-order. However, the type of transition is variant for 
pattern (a). It is found that for most of t' values the transion is of second-order, while when 
t' is larger than about 0.17 it becomes of first-order. The explanation for the first-order here 
is similar to that for the first-order suppression by i! at T = 0, which was presented in Ref. 
0- 

. 12 
0.1 
0.08 
^ 0.0 6 
0.04 
0.02 

0.05 0.1 0.15 0.2 
T 

FIG. 5. The optimal dimerization parameter 5* as a function of temperature T un- 
der various values of t' for pattern (a). The curves from right to left correspond to 
t' = 0, 0.12, 0.16, 0.168, 0.17, 0.171, 0.172, respectively. 
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FIG. 6. The same as Fig. |5| for pattern (b). The curves from right to left correspond to 
t' = 0, 0.1, 0.12, 0.15, 0.16, 0.168, 0.17, respectively. 
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To see the competition between the two dimerization patterns, one needs to compare 
the free energies F* for both patterns at their respective optimal dimerization values for 
each if and T. In the following we choose several t' values 0, 0.14, 0.168, 0.17 to show all 
different kinds of results, see Figs. [7|-|l0|. In each figure the optimal dimerization parameters 
and their corresponding free energies (per site) with change of temperature are compared 
for both patterns. We will discuss them in detail in the following. 

For the cases of if = 0, 0.14, 0.168, it is always found that the transition temperatures for 
both patterns are identical, i.e., = T p = T p , but the comparative results on free energies 
are different. For if — the free energy for pattern (b) F£ is always lower than that for 
pattern (a) F* in the whole temperature region T < T p (see Fig. [7]), which means that the 
dimerized state with pattern (b) is always prefered; while for if = 0.168 the result is exactly 
the reverse (see Fig. ||). The most interesting case is that for the intermediate t' = 0.14, 
see Fig. |8|. To more clearly compare the magnitudes of F* and F£, we have plotted their 
difference AF* = F£ — F* as a function of temperature in the lowest panel, which shows that 
_F 6 * is lower than F* in the intermediate temperature region and then becomes higher for low 
temperatures. Note that the absolute numerical error for the free energy in our calculations 
is less than 10 -8 , so the minor difference AF* is not an artifact. Therefore, in this case 
the transitions with decreasing temperature may be described as follows: at some critical 
temperature T c x the system goes through a Peierls instability into the dimerized state (b). 
When temperature continues to decrease down to another critical value T c 2 the system will 
turn from the dimerized state (b) into (a). That is to say, in the whole temperature region 
two successive Peierls transitions take place, as expected in the Introduction. 

In addition, for the case of t' = 0.17, the transition temperature T£ becomes different 
from Ti 1 and the former is larger, see Fig. [H]. In addition, the free energy F* is always lower 
than Thus the pattern (a) is always dominant once T < T?. 
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FIG. 7. The optimal values S* (upper panel) and the corresponding free energies F* (lower 
panel) as a function of T for t' = 0. The dashed lines are for pattern (a) and the solid lines are for 
pattern (b). 
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With the known results for several typical t' values, we now construct the phase diagram 
for different stable states in the full parameter space of t' and T. The result is presented 



in Fig. [IT], in which most of the above results have been summarized. From Fig. |TTJ the 
transition temperatures T£ (solid line) and (dotted line) with change of t', as well as 
the competition between two dimerization patterns are clearly seen. Note that the solid 
and dotted lines are identical for most values of t', and become separate only for a narrow 



region with large t' . The most exciting result seen from Fig. |TT] is that the double Peierls 
transitions, i.e., the system passing the states U— >D2— >D1 with decrease of temperature, 
may take place for t' in the region with approximately 0.137 <t'< 0.16. 




0.05 



0.2 



FIG. 11. The phase diagram for different stable states in the parameter space of t' and T. U, 
Dl, D2 represent the uniform state, dimerized state with pattern (a), dimerized state with pattern 
(b), respectively. The solid and dotted lines are the phase boundaries between U and Dl, and D2, 
respectively, when the two dimerization patterns are considered separately. Note that they are 
identical for most values of t'. The dashed line differentiates the stable state between Dl and D2. 
In the narrow area between the solid and dotted line only the Dl state survives. 



V. DISCUSSION AND CONCLUSION 

In the previous sections, extensive results about the Peierls instability, including the 
rich phase diagram in the parameter space of t' and T have been discovered. Now we 
discuss the possible implications of these results to real materials. Although our model 
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may be too simplified when fitting to real quasi-2D Peierls materials, the results obtained 
are still instructive for the understanding of some experimental findings. (Actually some of 
instructions have already been presented in the suitable context above.) For real materials, 
the n.n.n. hopping t' is often indispensable, and moreover, it can be changed by some 
external conditions, e.g., pressure. Depending on whether the ratio t'/t is enhanced or 
reduced by the pressure, the Peierls transition temperature is expected to decrease (and 
even vanish) or increase by it. Experimentally both of possibilities were observed in the 
real quasi-2D bronzes, see Ref. || and references therein. Our results also suggest that 
applying pressure (if leading to an enhancement of t'/t) may prevent the materials from the 
Peierls instability, and then possibly help them enter into the competing superconducting 
state. Interestingly for the bronze Li .9Mo 6 Oi7, it was found that the superconducting 
transition temperature is largely increased by pressure, associated with a sharp decrease of 
the transition temperature of a possible Peierls instability ||18|| . In particular, our results 



show two successive Peierls transitions with temperature, a phenomenon which was found 
in TlMo 6 Oi 7 an d (P0 2 )4(W0 3 )2 m Essentially, the theoretical and experimental 

results are similar, i.e., they both originate from the competition between two different 
lattice distortion patterns. Of course, the real situations are more complicated mainly 
because of intricate lattice distortions. Nevertheless, our work provides a good example in 
two dimensions to show double Peierls transitions at finite temperatures. Note that, within 
the intrinsic limitations (i.e., adiabatic phonons and limited lattice distortion patterns), all 
results presented here are exact in the thermodynamic limit. 

Coming back to the region of Fig. |TT| where successive transitions U— >D2-^D1 take place, 
we further discuss an interesting issue stated as follows. In this figure, it is obvious that the 
second transition D2— >D1 is of first-order. Are there some intermediate states which may 
smoothly connect D2 and Dl? Actually, as a more general dimerization pattern, the lattice 
distortions may be assumed in the way: 
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with a ratio r = u y /u x G [0, 1]. The patterns (a) and (b) correspond to the two limiting cases: 
r = 1 and 0, respectively. If this generalization is included in the calculation, preliminary 
results give the following scenario with decreasing temperature: first, the transition still 
happens between the uniform state and the dimerized state (b) (i.e., r = 0). This state 
persists for a range of T and then begins to turn continuously into other dimerized states 
with r > 0. Depending on the value of t', the r = 1 state, i.e., the dimerized state (a) may 
be quickly reached or may not be reached until T = 0. Now, the second transition between 
Dl and other dimerized states becomes of second-order. Detailed results need a more careful 
calculation. 

Finally, we briefly discuss the effects of electron correlations, e.g., the Hubbard on- 
site U at ground state. For t' = 0, it was found that the on-site Coulomb interaction is 
unfavorable to dimerization in two dimensions as soon as U is present ]10| , 20 1 . This is 



because the on-site U favors the appearance of antiferromagnetic (AF) order for the 2D 
half-filled model, while the dimerization tends to stabilize local spin singlets. Thus both, 
Peierls and AF instabilities, will compete with each other in the presence of finite r\ and U 
|20|| . On the other hand, the hopping t', which breaks the perfect nesting, will frustrate the 
Peierls instability as shown here and also the AF instability as shown by several authors 
19| , when each of them is discussed individually. Then it surely has nontrivial effects on 



the competition between these two instabilities when they are considered jointly. As far as 
the Peierls instability is concerned, it is argued that both t' and U may cancel in part their 
effects when acting simultaneously, although each of them separately tends to suppress it. 
This topic will be left for future investigation. 

In conclusion, the Peierls instabilities for a half-filled 2D t-t' model are thoroughly studied 
with consideration of two dimerization patterns (a) and (b). The effect of imperfect nesting 
introduced by t' on the Peierls instability, the properties of the dimerized ground state, as 
well as the competition between two dimerized states are completely discussed. We have 
found that the n.n.n. hopping t' will frustrate the Peierls instability for each of the dimerized 
states when considered separately. Moreover, when the two dimerized states are considered 



together, they will compete to result in a rich phase diagram in the parameter space of 
t' and T. Prominently, double Peierls transitions, that is, the system passing from the 
uniform state to one dimerized state and then to the other may take place with decrease of 
temperature. 
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